skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Search for: All records

Creators/Authors contains: "Johnson, James W"

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. Abstract Many nucleosynthetic channels create the elements, but two-parameter models characterized byαand Fe nonetheless predict stellar abundances in the Galactic disk to accuracies of 0.02–0.05 dex for most measured elements, near the level of current abundance uncertainties. It is difficult to make individual measurements more precise than this to investigate lower-amplitude nucleosynthetic effects, but population studies of mean abundance patterns can reveal more subtle abundance differences. Here, we look at the detailed abundances for 67,315 stars from the Apache Point Observatory Galactic Evolution Experiment (or APOGEE) Data Release 17, but in abundance residuals away from a best-fit two-parameter, data-driven nucleosynthetic model. We find that these residuals show complex structures with respect to age, guiding radius, and vertical action that are not random and are also not strongly correlated with sources of systematic error such as log ( g ) ,Teff, and radial velocity. The residual patterns, especially in Na, C+N, Mn, and Ce, trace kinematic structures in the Milky Way, such as the inner disk, thick disk, and flared outer disk. A principal component analysis suggests that most of the observed structure is low-dimensional and can be explained by a few eigenvectors. We find that some, but not all, of the effects in the low-αdisk can be explained by dilution with fresh gas, so that the abundance ratios resemble those of stars with higher metallicity. The patterns and maps we provide can be combined with accurate forward models of nucleosynthesis, star formation, and gas infall to provide a more detailed picture of star and element formation in different Milky Way components. 
    more » « less
    Free, publicly-accessible full text available April 25, 2026
  2. Abstract Metallicities of both gas and stars decline toward large radii in spiral galaxies, a trend known as the radial metallicity gradient. We quantify the evolution of the metallicity gradient in the Milky Way as traced by APOGEE red giants with age estimates from machine learning algorithms. Stars up to ages of ∼9 Gyr follow a similar relation between metallicity and Galactocentric radius. This constancy challenges current models of Galactic chemical evolution, which typically predict lower metallicities for older stellar populations. Our results favor anequilibrium scenario, in which the gas-phase gradient reaches a nearly constant normalization early in the disk lifetime. Using a fiducial choice of parameters, we demonstrate that one possible origin of this behavior is an outflow that more readily ejects gas from the interstellar medium (ISM) with increasing Galactocentric radius. A direct effect of the outflow is that baryons do not remain in the ISM for long, which causes the ratio of star formation to accretion, Σ ̇ / Σ ̇ in , to quickly become constant. This ratio is closely related to the local equilibrium metallicity, since its numerator and denominator set the rates of metal production by stars and hydrogen gained through accretion, respectively. Building in a merger event results in a perturbation that evolves back toward the equilibrium state on ∼Gyr timescales. Under the equilibrium scenario, the radial metallicity gradient is not a consequence of the inside-out growth of the disk but instead reflects a trend of declining Σ ̇ / Σ ̇ in with increasing Galactocentric radius. 
    more » « less
    Free, publicly-accessible full text available July 10, 2026
  3. ABSTRACT The scaling of the specific Type Ia supernova (SN Ia) rate with host galaxy stellar mass $$\dot{\text{N}}_\text{Ia} / \text{M}_\star \sim \text{M}_\star ^{-0.3}$$ as measured in ASAS-SN and DES strongly suggests that the number of SNe Ia produced by a stellar population depends inversely on its metallicity. We estimate the strength of the required metallicity dependence by combining the average star formation histories (SFHs) of galaxies as a function of their stellar mass with the mass–metallicity relation (MZR) for galaxies and common parametrizations for the SN Ia delay-time distribution. The differences in SFHs can account for only ∼30 per cent of the increase in the specific SN Ia rate between stellar masses of M⋆ = 1010 and 107.2 M⊙. We find that an additional metallicity dependence of approximately ∼Z−0.5 is required to explain the observed scaling. This scaling matches the metallicity dependence of the close binary fraction observed in APOGEE, suggesting that the enhanced SN Ia rate in low-mass galaxies can be explained by a combination of their more extended SFHs and a higher binary fraction due to their lower metallicities. Due to the shape of the MZR, only galaxies below M⋆ ≈ 3 × 109 M⊙ are significantly affected by the metallicity-dependent SN Ia rates. The $$\dot{\text{N}}_\text{Ia} / \text{M}_\star \sim \text{M}_\star ^{-0.3}$$ scaling becomes shallower with increasing redshift, dropping by factor of ∼2 at 107.2 M⊙ between z = 0 and 1 with our ∼Z−0.5 scaling. With metallicity-independent rates, this decrease is a factor of ∼3. We discuss the implications of metallicity-dependent SN Ia rates for one-zone models of galactic chemical evolution. 
    more » « less
  4. ABSTRACT We model the stellar abundances and ages of two disrupted dwarf galaxies in the Milky Way stellar halo: Gaia-Sausage Enceladus (GSE) and Wukong/LMS-1. Using a statistically robust likelihood function, we fit one-zone models of galactic chemical evolution with exponential infall histories to both systems, deriving e-folding time-scales of τin = 1.01 ± 0.13 Gyr for GSE and $$\tau _\text{in} = 3.08^{+3.19}_{-1.16}$$ Gyr for Wukong/LMS-1. GSE formed stars for $$\tau _\text{tot} = 5.40^{+0.32}_{-0.31}$$ Gyr, sustaining star formation for ∼1.5–2 Gyr after its first infall into the Milky Way ∼10 Gyr ago. Our fit suggests that star formation lasted for $$\tau _\text{tot} = 3.36^{+0.55}_{-0.47}$$ Gyr in Wukong/LMS-1, though our sample does not contain any age measurements. The differences in evolutionary parameters between the two are qualitatively consistent with trends with stellar mass M⋆ predicted by simulations and semi-analytic models of galaxy formation. Our inferred values of the outflow mass-loading factor reasonably match $$\eta \propto M_\star ^{-1/3}$$ as predicted by galactic wind models. Our fitting method is based only on Poisson sampling from an evolutionary track and requires no binning of the data. We demonstrate its accuracy by testing against mock data, showing that it accurately recovers the input model across a broad range of sample sizes (20 ≤ N ≤ 2000) and measurement uncertainties (0.01 ≤ σ[α/Fe], σ[Fe/H] ≤ 0.5; $$0.02 \le \sigma _{\log _{10}(\text{age})} \le 1$$). Due to the generic nature of our derivation, this likelihood function should be applicable to one-zone models of any parametrization and easily extensible to other astrophysical models which predict tracks in some observed space. 
    more » « less
  5. ABSTRACT We derive empirical constraints on the nucleosynthetic yields of nitrogen by incorporating N enrichment into our previously developed and empirically tuned multizone galactic chemical evolution model. We adopt a metallicity-independent (‘primary’) N yield from massive stars and a metallicity-dependent (‘secondary’) N yield from AGB stars. In our model, galactic radial zones do not evolve along the observed [N/O]–[O/H] relation, but first increase in [O/H] at roughly constant [N/O], then move upward in [N/O] via secondary N production. By t ≈ 5 Gyr, the model approaches an equilibrium [N/O]–[O/H] relation, which traces the radial oxygen gradient. Reproducing the [N/O]–[O/H] trend observed in extragalactic systems constrains the ratio of IMF-averaged N yields to the IMF-averaged O yield of core-collapse supernovae. We find good agreement if we adopt $$y_\text{N}^\text{CC}/y_\text{O}^\text{CC}=0.024$$ and $$y_\text{N}^\text{AGB}/y_\text{O}^\text{CC} = 0.062(Z/Z_\odot)$$. For the theoretical AGB yields we consider, simple stellar populations release half their N after only ∼250 Myr. Our model reproduces the [N/O]–[O/H] relation found for Milky Way stars in the APOGEE survey, and it reproduces (though imperfectly) the trends of stellar [N/O] with age and [O/Fe]. The metallicity-dependent yield plays the dominant role in shaping the gas-phase [N/O]–[O/H] relation, but the AGB time-delay is required to match the stellar age and [O/Fe] trends. If we add ∼40 per cent oscillations to the star formation rate, the model reproduces the scatter in the gas phase [N/O]–[O/H] relation observed in external galaxies by MaNGA. We discuss implications of our results for theoretical models of N production by massive stars and AGB stars. 
    more » « less
  6. Abstract We investigate the [X/Mg] abundances of 16 elements for 82,910 Galactic disk stars from GALAH+ DR3. We fit the median trends of low-Ia and high-Ia populations with a two-process model, which describes stellar abundances in terms of a prompt core-collapse and delayed Type-Ia supernova component. For each sample star, we fit the amplitudes of these two components and compute the residual Δ[X/H] abundances from this two-parameter fit. We find rms residuals ≲0.07 dex for well-measured elements and correlated residuals among some elements (such as Ba, Y, and Zn) that indicate common enrichment sources. From a detailed investigation of stars with large residuals, we infer that roughly 40% of the large deviations are physical and 60% are caused by problematic data such as unflagged binarity, poor wavelength solutions, and poor telluric subtraction. As one example of a population with distinctive abundance patterns, we identify 15 stars that have 0.3–0.6 dex enhancements of Na but normal abundances of other elements from O to Ni and positive average residuals of Cu, Zn, Y, and Ba. We measure the median elemental residuals of 14 open clusters, finding systematic ∼0.1–0.4 dex enhancements of O, Ca, K, Y, and Ba and ∼0.2 dex depletion of Cu in young clusters. Finally, we present a restricted three-process model where we add an asymptotic giant branch star (AGB) component to better fit Ba and Y. With the addition of the third process, we identify a population of stars, preferentially young, that have much higher AGB enrichment than expected from their SNIa enrichment. 
    more » « less
  7. Abstract We measure abundances of 12 elements (Na, Mg, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni) in a sample of 86 metal-poor (−2 ≲ [Fe/H] ≲ −1) subgiant stars in the solar neighborhood. Abundances are derived from high-resolution spectra taken with the Potsdam Echelle Polarimetric and Spectroscopic Instrument on the Large Binocular Telescope, modeled using iSpec and MOOG. By carefully quantifying the impact of photon-noise (<0.05 dex for all elements), we robustly measure theintrinsicscatter of abundance ratios. At fixed [Fe/H], the rms intrinsic scatter in [X/Fe] ranges from 0.04 (Cr) to 0.16 dex (Na), with a median of 0.08 dex. Scatter in [X/Mg] is similar, and accounting for [α/Fe] only reduces the overall scatter moderately. We consider several possible origins of the intrinsic scatter with particular attention to fluctuations in the relative enrichment by core-collapse supernovae (CCSN) and Type Ia supernovae and stochastic sampling of the CCSN progenitor mass distribution. The stochastic sampling scenario provides a good quantitative explanation of our data if the effective number of CCSN contributing to the enrichment of a typical sample star isN∼ 50. At the median metallicity of our sample, this interpretation implies that the CCSN ejecta are mixed over a gas mass ∼6 × 104Mbefore forming stars. The scatter of elemental abundance ratios is a powerful diagnostic test for simulations of star formation, feedback, and gas mixing in the early phases of the Galaxy. 
    more » « less
  8. Abstract We report the first direct measurement of the helium isotope ratio, 3 He/ 4 He, outside of the Local Interstellar Cloud, as part of science-verification observations with the upgraded CRyogenic InfraRed Echelle Spectrograph. Our determination of 3 He/ 4 He is based on metastable He i * absorption along the line of sight toward Θ 2 A Ori in the Orion Nebula. We measure a value 3 He/ 4 He = (1.77 ± 0.13) × 10 −4 , which is just ∼40% above the primordial relative abundance of these isotopes, assuming the Standard Model of particle physics and cosmology, ( 3 He/ 4 He) p = (1.257 ± 0.017) × 10 −4 . We calculate a suite of galactic chemical evolution simulations to study the Galactic build up of these isotopes, using the yields from Limongi & Chieffi for stars in the mass range M = 8–100 M ⊙ and Lagarde et al. for M = 0.8–8 M ⊙ . We find that these simulations simultaneously reproduce the Orion and protosolar 3 He/ 4 He values if the calculations are initialized with a primordial ratio 3 He / 4 He p = ( 1.043 ± 0.089 ) × 10 − 4 . Even though the quoted error does not include the model uncertainty, this determination agrees with the Standard Model value to within ∼2 σ . We also use the present-day Galactic abundance of deuterium (D/H), helium (He/H), and 3 He/ 4 He to infer an empirical limit on the primordial 3 He abundance, 3 He / H p ≤ ( 1.09 ± 0.18 ) × 10 − 5 , which also agrees with the Standard Model value. We point out that it is becoming increasingly difficult to explain the discrepant primordial 7 Li/H abundance with nonstandard physics, without breaking the remarkable simultaneous agreement of three primordial element ratios (D/H, 4 He/H, and 3 He/ 4 He) with the Standard Model values. 
    more » « less
  9. ABSTRACT We test the hypothesis that the observed first-peak (Sr, Y, Zr) and second-peak (Ba) s-process elemental abundances in low-metallicity Milky Way stars, and the abundances of the elements Mo and Ru, can be explained by a pervasive r-process contribution originating in neutrino-driven winds from highly magnetic and rapidly rotating proto-neutron stars (proto-NSs). We construct chemical evolution models that incorporate recent calculations of proto-NS yields in addition to contributions from asymptotic giant branch stars, Type Ia supernovae, and two alternative sets of yields for massive star winds and core-collapse supernovae. For non-rotating massive star yields from either set, models without proto-NS winds underpredict the observed s-process peak abundances by 0.3–$$1\, \text{dex}$$ at low metallicity, and they severely underpredict Mo and Ru at all metallicities. Models incorporating wind yields from proto-NSs with spin periods P ∼ 2–$$5\, \text{ms}$$ fit the observed trends for all these elements well. Alternatively, models omitting proto-NS winds but adopting yields of rapidly rotating massive stars, with vrot between 150 and $$300\, \text{km}\, \text{s}^{-1}$$, can explain the observed abundance levels reasonably well for [Fe/H] < −2. These models overpredict [Sr/Fe] and [Mo/Fe] at higher metallicities, but with a tuned dependence of vrot on stellar metallicity they might achieve an acceptable fit at all [Fe/H]. If many proto-NSs are born with strong magnetic fields and short spin periods, then their neutrino-driven winds provide a natural source for Sr, Y, Zr, Mo, Ru, and Ba in low-metallicity stellar populations. Conversely, spherical winds from unmagnetized proto-NSs overproduce the observed Sr, Y, and Zr abundances by a large factor. 
    more » « less